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Systems containing few Fermions (e.g., electrons) are of 
great current interest. Fluorescence occurs when electrons 
drop from one level to another without changing spin. Only 
electron gases in a state of equilibrium are considered. When 
the system may exchange electrons with a large reservoir, 
the electron-gas fluorescence is easily obtained from the well- 
known Fermi-Dirac distribution. But this is not so when the 
number of electrons in the system is prevented from vary- 
ing, as is the case for isolated systems and for systems that 
are in thermal contact with electrical insulators such as dia- 
mond. Our accurate expressions rest on the assumption that 
single-electron energy levels are evenly spaced, and that en- 
ergy coupling and spin coupling between electrons are small. 
These assumptions are shown to be realistic for many systems. 
Fluorescence from short, nearly isolated, quantum wires is 
predicted to drop abruptly in the visible, a result not pre- 
dicted by the Fermi-Dirac distribution. Our exact formulas 
are based on restricted and unrestricted partitions of integers. 
The method is considerably simpler than the ones proposed 
earlier, which are based on second quantization and contour 
integration. 
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I. INTRODUCTION 

A number of remarkable experiments involving few 
electrons in semiconductors and free-space, metal parti- 
cles, and spin-1/2 atoms at low temperatures, have been 
recently reported [0, Only electrons are considered 
below. These collections of electrons may be isolated or 
in thermal contact with the environment, but, in any 
event, the number of particles is constant. The Fermi- 
Dirac (FD) distribution holds when electrons may be 
freely exchanged with a large reservoir (grand-canonical 
ensemble), but is inaccurate for the systems considered. 
The present paper provides simple and accurate formulas 
for electron occupancy and fluorescence for evenly-spaced 
single-electron energy levels. Spontaneous emission is 
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supposed to be weak enough not to perturb importantly 
the system state of equilibrium Only vanishingly 
small Coulomb interaction between electrons is con- 
sidered. The time required for the system to reach equi- 
librium is not needed because averagings are performed 
over unlimited time scales. Quantum Optics effects, such 
as resonance fluorescence or super-radiance, will not be 
considered. 

The assumption of evenly-spaced single-electron en- 
ergy levels is not as restrictive as one may think at 
first. Consider indeed one-dimensional devices such as 
the quantum wires employed in modern laser diodes [Q. 
If the wire is uniform over its length and the valence and 
conduction bands are parabolic in the energy-momentum 
space, the energy spacing, e, between adjacent levels is 
not a constant. However, the variations of e may be ne- 
glected near the Fermi level as long as the temperature 
is not too high. For zero-band-gap semiconductors such 
as Pbo.84Sno.i6Te, the energy-momentum relationship is 
linear rather than parabolic, and e is exactly a constant. 
Level spacings in small irregular metal particles (with 
a size on the order of 10 nm) are nearly uniform as a 
consequence of the mechanism of level "repulsion" . The 
probability that adjacent levels be separated by e is, for 
the appropriate ensemble, of the form: e*exp(— e^), a 
sharply peaked function of e [Q. As is well-known, the 
Landau levels that describe electron motion in uniform 
magnetic fields Q are evenly spaced. These levels are 
highly degenerate, but the coupling between degenerate 
states (expressing the drift of electrons through the mag- 
netic field lines) may be neglected over some period of 
time. Likewise, two or three-dimensional harmonic oscil- 
lators, modeling for example the confinement of electrons 
in traps, exhibit degenerate evenly-spaced levels. Our 
approach may be generalized to degenerate levels. It is 
appropriate to mention also that the density of states of 
(two-dimensional) quantum wells with parabolic bands 
is a constant within a sub-band. This implies that the 
energy-level spacing is constant on the average, though 
not exactly. 

The amount of light spontaneously emitted by elec- 
tronic systems depends on the optical-mode density of 
state, which is different for free-space, low-dimensional 
structures, or photonic band gap materials. It is not the 
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purpose of the present paper to discuss such problems. 
Because aU comparisons are made at the same optical 
wavelength, terms depending on wavelength only (essen- 
tially the optical-mode density of state) may be dropped. 
The quantity that we calculate is the probability that the 
system exhibits an electron at level k and a hole (no elec- 
tron) at level k' , with k — k' = fku / e = d, where h denotes 
the Planck constant divided by 27r and lo the angular op- 
tical frequency of observation. In pure semiconductors, 
electron-momentum conservation entails that transitions 
may occur only symmetrically with respect to the Fermi 
level, implying that k = {l + d)/2, k' = (1 - d)/2, with 
odd d, if k — labels the zero-temperature top electron. 
We will particularly consider the case where this condi- 
tion holds. 

Let us recall that, in Statistical Mechanics, isolated 
systems are described by microcanonical ensembles, sys- 
tems that may exchange energy but not particles with a 
reservoir are described by canonical ensembles, and sys- 
tems that may exchange both energy and particles with 
a reservoir are described by grand-canonical ensembles. 
The Fermi-Dirac (FD) distribution Q is applicable to fi- 
nite systems only in the latter case. That is, the FD dis- 
tribution is invalid for isolated systems and for systems 
that are in contact with an electrical insulator such as 
diamond. In the present context, "finite" (or "small", or 
"mesoscopic" ) means that fc^T, where T denotes the ab- 
solute electron gas temperature, is not necessarily large 
compared with the adjacent-level energy spacings e. In 
grand canonical ensembles, the fluorescence is propor- 
tional to the product of the probabilities that the upper 
level be occupied and that the lower level be empty. In 
canonical ensembles, it turns out that fluorescence is pro- 
portional to the difference between lower and upper levels 
occupancies. This is apparently a new result. 

Exact formulas for level occupancy in finite single-spin 
systems with evenly-spaced levels in contact with a heat 
bath have been reported before §, §, §, g, §. Our 
method ||] consists of first enumerating the microstates 
of isolated systems. Subsequent averaging provides ex- 
pressions for the canonical occupancies. This method is 
considerably simpler than the second-quantization meth- 
ods and integral-transformation formulas employed in [Q , 
The present paper generalizes the results reported in 
p I to account for the fact that some electrons may change 
spin in the course of time. Simple, exact formulas for flu- 
orescence are obtained for isolated systems and canonical 
ensembles. For arbitrary level energies and the canonical 
ensemble, known recurrence formulas are satisfied by 
our more special, but explicit, forms. 

The FD distribution is derived in Section 2 to set up 
the notation and to recall why, in grand-canonical ensem- 
bles, spin is properly accounted for by a two-fold degen- 
eracy. This method of accounting for the electron spin, 
however, is invalid in microcanonical or canonical ensem- 
bles because unpaired electrons may be in either one of 
two distinguishable states. Finite systems that exchange 
energy but not particles with a reservoir are considered 



in Section 3 and 4. In Section 3, the electron spin is 
supposed to be strictly preserved and, for simplicity, it 
is supposed that there are as many electrons with spin 
up and with spin down. Formulas valid for single-spin 
electrons need only be multiplied by factors of 2 in that 
situation. In Section 4, electrons arc allowed to change 
spin in the course of time (but not during a spontaneous 
emission event). For the sake of clarity, only essential 
formulas relating to canonical ensembles are given in the 
main text, detailed derivations and intermediate results 
being relegated to Appendices A and B. 

The purpose of Appendix |^ is to explain why the total 
number W of microstates of isolated systems is related to 
the partitions of integers, for single-spin-state electrons. 
It is shown that the number m{k) of microstates whose 
level k is occupied is simply related to W. The num- 
ber m(k; k') of microstates whose level k is occupied and 
level k' is empty is shown to be simply related to m{k), 
and thus to W. Averaging, with the Boltzmann factor 
as a weight, provides the corresponding formulas in the 
canonical ensemble. Appendix ^ explains how the pos- 
sibility that electrons may change spin in the course of 
time is accounted for. 



II. THE FERMI-DIRAC DISTRIBUTION 

The probability p that a system in thermal and elec- 
trical contact with a large medium contains elec- 
trons with spin up, N~ electrons with spin down, and 
energy U, is proportional to the corresponding number 
of medium states (subscripts m) conveniently written as 
exp (S'm(A'+, iV^, [/„)) with Sm the system entropy. If 
the system-medium contact is very weak, energies as well 
as particle numbers add up. A first-order expansion of 
Sm with respect to its arguments then gives 

p{N+,N-,U)^Ceii-p{~-aN+ -aN- -f3U) (1) 

where C denotes a constant and 
95^, 
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Here, /3 = l/ksT, where T denotes the temperature, and 
fi the Fermi level. A single /x value occurs because the 
medium behavior is the same for electrons of opposite 
spins. 

It follows from Eq. (|l|) that the probabilities p^*'-' that a 
nondegenerate state of energy efe is unoccupied, p^^^ that 
it is occupied by an electron of either spin, and p^-^^^ that 
it is occupied by two electrons, are respectively 
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with z = exp(— a — /3efe). Normalization [jP*^ + 2p(^^ 



= 1) gives C = 1/(1 + zY . The occupancy (average 
number of electrons) npD = 2p(^) + ^^P'^ of level k is 
therefore 



z-i + 1 



exp(/3(efc - ^)) + 1 
2 



92 + 1 



(4) 



where g = exp(— /3) denotes the Boltzmann factor. In 
the last expression it is assumed that ek — k, where k 
denotes any relative integer. The separation e between 
adjacent-level energies is taken as the energy unit, with a 
typical value for 1 /xm-long quantum wires e = 1 meV. In 
the last expression in Eq. (^, fc = 1 labels the level just 
above the top electron at T = K, and we have set /i = i. 

Note that for large fc-values, npoiq] k) « 2q^^^/'^. 

The average system energy added on top of the T — 
K energy is obtained by summing up the occupancy 
over all levels, and subtracting a similar sum for the T ~ 
K distribution. The result is lol 



3=1,3.... ^ 



(5) 



In grand canonical ensembles, occupancies at different 
levels are independent. For a single spin state, this means 
that the probability that level k is occupied and level k' 
is empty is the product of level k occupancy and [1— level 
k' occupancy]. When the two spin-states are considered, 
we obtain 



LFD{q;k,k') = 2 



npoiq, k) 



1 - 



nFD{q,k') 



(6) 



Fluorescence may indeed occur for 8 out of the 16 pos- 
sibilities of occupancy of levels k and k' (no electron, 
spin-up electron, spin-down electron, or two electrons for 
each of the two levels). Because occupancies are inde- 
pendent, the sum of the probabilities that fluorescence 
events occur is indeed found to be given by Eq. (^). 

If the electron-momentum conservation law is enforced, 
we have k = {1 + d)/2, fc' = (1 - d)/2, and the FD 
fluorescence in Eq. m) reads after rearranging 



Lpniq^d) 



{q-d/^ + iy 



(7) 



Thus, the fluorescence in the grand-canonical ensemble 
is given by a simple function of temperature T and angu- 
lar optical frequency of observation lo. We will see that 
the canonical ensemble fluorescence is given by a simple 
series. 



III. FLUORESCENCE WITHOUT SPIN FLIP 

For single-spin electrons, the occupancy in isolated sys- 
tems is given by a simple formula reported in Jp|. The 
proof, omitted in i s giv en in Appendix ^oi the 
present paper (see Eq. (A_12)). If the energy added to 
the system is denoted by r, the number W{r) of configu- 
rations of the system is equal here to the number p{r) of 
partitions of r. Indeed, microstates may be obtained by 
shifting electrons upward from their T = K locations 
by non-increasing steps that sum up to r. Let us recall 
that a partition of r is a set of positive integers summing 
up to r. For example (2,1,1) is a partition of 4. The 
number p(4) of partitions of 4 equals 5. By convention, 
p(0) = 1 and p{r) = if r < 0. 

To illustrate the difference existing between the exact 
result and the Fermi-Dirac distribution, let us note that, 
for any microstate, the energy separation between the top 
electron and the lowest hole cannot exceed re, where e 
denotes as before the adjacent-level energy spacing. Ac- 
cordingly, the fluorescence drops abruptly at an angular 
optical frequency lj — re/h. When r 3> 1 , a system tem- 
perature T may be defined ksT e\/6r/7r. As an 
example, room-temperature isolated systems with e = 1 
meV should not emit visible light according to the exact 
formula, while some faint light is expected according to 
the FD distribution. 

In the present section, it is supposed that the electron 
spins are preserved in the course of time, and that there 
are as many electrons with spin up and spin down, for 
simplicity. It then suffices to multiply the expressions 
for single-spin electron occupancy, average energy, and 
fluorescence, given in Appendix H, by factors of 2. 

The occupancy reads, according to Eq. (A13) 



n„(g;fc) = -2 ^ ^^ly q3k+oU-i)/2 (g) 
3=1,2... 

Note that, for large fc-values, nu{q;k) w 2q'^ so that 
n.a{q; k)/nFD{q\ k) — if the expression for npoiq] k) 
in Eq. (^ is used. 

The average added energy reads m 



2 E 

j=l,2.. 



- 1 



(9) 



We first compare in Fig. |^a the occupancy in isolated 
systems with the FD occupancy. The forme r is obtained 
by multiplying ni(r;k) given in Eq. (A12) by a factor 
of 2 to account for the two spin states. The FD occu- 
pancy is given in Eq. (^, with q expressed in terms of 
the average energy r with the help of Eq. (||). Note that, 
below 0.1, isolated-system occupancies are .smaller than 
FD occupancies. 

Consider next the case where the system is in contact 
with a heat bath. The ratio of the exact occupancy in Eq. 
(^ (where q is expressed in terms of the average energy 
r with the help of Eq. (||)) and the FD occupancy (Eqs. 
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®' (H))' represented in Fig. |i|b as a function of the 
FD occupancy for various values of the average energy 
r (namely, r — 6, 60, and 600). Below 0.1, the exact 
occupancy exceeds the FD occupancy. Figure ^ shows 
that, when the comparison is made at equal temperatures 
T (instead of equal average energies), the opposite occurs. 
In Fig. H, we have chosen to represent the occupancy 
ratios as functions of the level number k (referred to the 
Fermi level) instead of the FD occupancy. 

To evaluate fluorescence, we need to know the number 
77i(r; fc; k') of microstates of added energy r having an 
electron at level k and none at level k' = k — d, where 
d = hw/t. Appendix ^ shows that this quantity is easily 
expressible in terms of the numbers m(r; k) of microstates 
having an electron at level k (irrespectively of other state 
occupancies). Averaging the result over r with as a 
weight to account for energy fluctuations in the presence 
of a heat bath, the fluorescence is found to be (see Eq. 



Lu{q;k,k') 



nu{q\k') - nuiq;k) 



q-'^ - 1 



(10) 



where ri„(g; k) is given in Eq. 

If the law of electron-momentum conservation is en- 
forced, the fluorescence reads 



Luiq;d) = i„(q; 



l+d l-d 



(11) 



Recall that q = exp(— e/ZcsT) (where e is typically 
1 meV and room temperature corresponds to fc^T = 
26 meV), and d = huj/e. The fluorescence ratio: 
Lu{q;d)/LpD{q;d) according to Eq. (0) and Eq. (11), 
is represented in Fig. M function of d for different 

temperatures. It is interesting that canonical and grand- 
canonical fluorescences almost coincide at small wave- 
lengths even though the occupancies are quite different 
in that limit. 



IV. FLUORESCENCE WITH SPIN FLIP 

Electron spins are now allowed to vary in the course of 
time (but not during a fluorescence event) . The numbers 
iV+ and N~ of electrons with spin up and spin down, re- 
spectively, may fluctuate, but their sum + = N 
remains constant if the system is isolated or in contact 
with an electrically-insulating heat sink such as diamond. 
The occupancy and fluorescence for coupled-spin elec- 
trons is derived from previous expressions through a suc- 
cession of partitionings and averagings. Because the de- 
tails are lengthy, they are relegated to Appendix ^. Re- 
markably, many summations can be performed in closed 
form so that the final result is simple. 

The occupancy reads 



■ °° „ '7"' n^i (q\ k + i) 



+ 00 



(12) 



where nu{q;k) is given in Eq. (||). Comparisons with 
the FD distributions are exemplified in Figs. |^d and ||b. 
With the help of theta functions |lj] the average added 
energy may be written as a simple sum (see Eq. ( [B5|) ) 



rc{q) = 2 



E 

3 = 1.2.. 



q-1 - 1 



(13) 



The fluorescence reads (see Eqs. (|A14|) and (|B6|)) 



Lc{q] d) 



+ 00 



(14) 



where Lu{q;k,k') is given in Eq. (|T^. Fluorescence is 
illustrated in Fig. Note that the exact result is closer 
to the FD result when electrons are allowed to change 
spin in the course of time. 



V. CONCLUSION 

When a (possibly small) system is in thermal and 
electrical contact with a large medium such as a piece 
of copper, the average number of electrons occupying 
some energy level (occupancy) is twice the value given by 
Fermi-Dirac (FD) formula. The fluorescence (light spon- 
taneously emitted without electron-spin flip), defined in 
terms of the probability that an electron at level k may 
drop to level k — fuv/e (w denotes the angular optical 
frequency of observation), is also easily obtained. But 
when the system is isolated, or in thermal contact with 
an electrical insulator, electron occupancies are given by 
different expressions. Because modern electronics often 
employ short quantum wires supported by diamond heat 
sinks, it is important to have at our disposal precise ex- 
pressions for occupancy and fluorescence in such situa- 
tions. The expressions obtained in this paper were illus- 
trated by comparison with the FD results. We considered 
the case where the electron spins are strictly maintained 
in the course of time (Section 3), and the case where spin- 
flip is allowed (Section 4). We found, for example, that 
small FD occupancies should be multiplied by approxi- 
mately exp(-e/2fcBT) « 0.22 if e = 1 meV and T = 4 K, 
a factor that differs significantly from unity. But, unex- 
pectedly, the fluorescence turns out to be given rather 
accurately by the FD distribution. 

Our mathematical approach is based on a direct enu- 
meration of the microstates, and the results are expressed 
in terms of the number of partitions of integers. This 
method is considerably simpler than those previously re- 
ported for similar models, both conceptually and alge- 
braically. A computer simulation has given results that 
are in very good agreement with the analytical formulas 
reported in this paper. 

It is our intention to report in the future analytical 
and numerical results relating to mesoscopic laser-diode 
light fluctuations. A preliminary step consists of consid- 
ering single-mode cavities incorporating the electron gas 
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at thermal equilibrium, with one or two bands of states 
(for a single band, see Appendix ^ of the present paper) . 
The intraband Auger effect and the stimulated transi- 
tions may be introduced at that stage. Next, the proba- 
bility that low-lying electrons be promoted to high-lying 
levels by the action of a (quiet or fluctuating) pump, and 
the probability that light quanta be absorbed, are intro- 
duced. At low power, our simulation gives output light 
fluctuations that agree very well with elementary laser- 
noise theory predictions ||l5[. At high power, new effects 
(temperature fluctuations, spectral-hole burning, statis- 
tical fluctuations of the optical gain . . .) occur that are 
difhcult to handle analytically ||T^. The analytical for- 
mulas reported in this paper are helpfull to assess the 
accuracy of the simulation in special situations. 



APPENDIX A: OCCUPANCY AND 
FLUORESCENCE FOR SINGLE-SPIN-STATE 
ELECTRONS. ARBITRARY ENERGY LEVELS 

We are only concerned in the main text with one- 
electron energy levels {efc} = Z. Rigorous occupancy 
formulas are obtained in the present Appendix by con- 
sidering first arbitrary {efc}. Eventually the number of 
levels is allowed to go to infinity. 

Consider an isolated system whose nondegenerate one- 
electron level energies are, in increasing order, ei, £2, 
. . . efc, . . . es, with N < B single-spin electrons. Accord- 
ing to the Pauli principle each level may be occupied by 
only or 1 electron. The system energy U is therefore 
the sum of N of the efc . The purpose of this Appendix is 
to evaluate: 

1. The number W{N, U) of possible ways of obtaining 
some given U (number of microstates, or "statisti- 
cal weight" ) . 

2. The number ■m{N, U; k) of microstates whose level 
k is occupied. The occupancy n{N, U;k) of level k 
is defined as m{N, U; k)/W{N, U). 

3. The number m{N, U ; fc; k') of microstates whose 
level k is occupied and level k' is empty. The 
fluorescence L{N, U; k; k') emitted by electrons 

level k' is defined as 



dropping from level k to 
m{N,U-k;k')/W{N, U). 



These evaluations will be presented in reversed or- 
der. Let us first relate the number m{N, U ; fc; k') of mi- 
crostates whose level k is occupied and level k' is empty, 
to the numbers m{N, U; k) defined above. For each mi- 
crostate, let the electron at level k be transferred to the 
lower empty level k' . The number N of electrons is unaf- 
fected but the total energy gets reduced from U to U — d 
where d = efc — efc' > 0, and the roles of k and k' are 
reversed. The following equality 



therefore holds. Now notice that, 

m{N, U\ k) - m{N, U; k; k') -- 
m{N, [/; k') - m{N, U; k'; k) 



(A2) 



because the two sides of the above equation count mi- 
crostates whose levels k an d k' are both occupied. W hen 
the expression in Eq. (A2) is introduced in Eq. ([Al[ ) it- 
eration gives a recurrence relation for m{N, U ; k\ k') that 
reads 



m{N, U; k; k') = 

{m{N,U--jd;k') 



m{N,U - jd;k)) (A3) 



j = l,2... 



The above series terminates when the total energy van- 
ishes, i.e., when jd exceeds U. 

Consider next the m{N, U; k) microstates whose levels 
k are occupied, and remove these electrons. The same 
number of microstates is obtained, with A'^ — 1 electrons, 
total energy U — ek, and no electron at level k. The 
number of these new microstates may be written as the 
difference between the total number of microstates and 
the number of microstates whose level k is occupied. We 
have therefore the identity 



m{N, U; k) = W{N - 1, [/ - efc) - m{N - 1, C/ - 



efc; k) 
(A4) 



After a sufficient number of iterations, either the energy 
or the number of electrons becomes negative and the last 
term vanishes. The quantity m{N, U; k) may therefore 
be written as a finite sum 



m{N,U;k) 



J2 {-iyW{N-j,U-jek) (A5) 



i=i,2.. 



The series terminates when j exceeds either N or U/ek- 
We have the obvious identity (the number of occupied 
states for the whole set of microstates being written in 
two different manners) 



NWiN, U)=Y^ m{N, U; k) 



(A6) 



k>l 



A system in contact with a heat bath at tempera- 
ture T is described by the canonical ensemble. Let 
us define as in the main text q = exp(— /3), where 
P = l/ksT. The so-called partition function Z{N,q) 
is the sum over U of q^W{N, U), and the average energy 
is iq/Z)dZ{N,q)/dq. When both sides of Eq. (|A|) are 



multiplied by q and summed over U we obtain, using 
Eq. 



miN, U; k- k') = m{N, U - d; k'; k) 



(Al) 



fe=l,2.. 
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Y^qU-i^>'W{N-3,U~jek) 
u 

1 ^ {-iyZ{l,q^)Z{N~j,q) (A7) 



j=l,2...N 

Indeed, for a single electron {N — 1), U may only take 
one of the values and the statistical weight W is unity. 
It follows that Z{l,q) is the sum over k f rom 1 to i? of 
q*^*. Note that Z{0,q) = 1. Equation (A7) was obtained 
earlier from a less direct proof. 

The occupancy (average number of electrons) n(q; k) 
of level k, is equal to the sum over U of q^mjN, U- k), 
divided by Z{N, q), where Tn{N, U;k) is given in Eq. ( A5) 
and Z(N,q) is defined in Eq. (A7) from a recurrence 
relation. Thus 



n{q;k) 



-Ej^i,2...i-^y<l''''Z{N-j,q) 
Z{N,q) 



(A8) 



Expression (|A8|) was reported before [ |10| . 

The probability that level k be occupied and that level 
k' be empty at temperature T is similarly obtained by 
summing q^m{N, U; k; k') over [/, and dividing the result 
by Z{N,q), where m{N,U]k;k') is given in Eq. ([A^). 
The result of the summation may be expressed as the 
difference of the lower and upper occupancies, according 
to 



L{q;k,k') 



n{q;k') -n{q;k) 
q^k'-^k _ 1 



(A9) 



where the occupancy n{q;k) is given in Eq. (A8). Thus 
the fluorescence is equal to the difference between the 
occupancies at k' and k multiplied by the Bose function. 

Let us now specialize the above formulas for the case 
where = k, k = 1 . . . B. Considering the displace- 
ment of the N electrons from their least-energy locations 
(fc = 1 to N) beginning to the one on top, we observe that 
W{N,U) is the number p{P;N,r) of partitions of the 
added energy r = U — N{N + l)/2 into at most N parts, 
none of them exceeding P = B — N . Note that the num- 
bers p{P] N, r) may be obtained from a recurrence rela- 
tion ^,p{P;N,r)-p{P;N~l,r) ^ p{P -l;N,r- N). 
Let us now change slightly our notation, letting /c = 
denote the top electron in the least energy configuration, 
and let us employ the added energy r instead of the total 
energy U as an argument. Equation (|A5D reads 



(AlO) 



m{N,r-k) = - {-ly 

= 1.2... 

xp{P + j;N- j, r-jk- j{j - l)/2) 



If r does not exceed N and P, it is intuitive that 
p{P; N, r) — p{r), where p{r) denote s the number of un- 
restricted partitions of r. Equation ( |A10| ) then simplifies 
to 

m(r;fc) = - ^ {-iyp{r~jk~jU-l)/2) (All) 

3 = 1,2... 



This expression was reported (for the first time to our 
knowledge), in |6|. If and P are infinite ({e/c} = Z), 
Eq. ( All] ) holds for any finite value of r and the corre- 
sponding singie-spin-state occupancy of an isolated sys- 
tem is 



ni{r; k) = m{r; k) /p{r) 



(A12) 



Averaging the numerator and denominator of above ex- 
pression with q^ as a weight with r from to cx), gives 
the canonical occupancy H] 



n{q;k) 



E 

i=i,2.. 



(-1) 



lYqjk+jij- 



l)/2 



(A13) 



Finally the expression in Eq. (A9) simplifies in the 
present situation to 



L{q;k,k') ^ 



n(q; k') - n{q; k) 

„k'~k _ I 



(A14) 



We have set up a Monte Carlo simulation program that 
enables us to recover previous analytical expressions. For 
the case of isolated systems, a constant probability per 
unit time is ascribed to level-changing events that pre- 
serve energy. The system eventually reaches a state of 
equilibrium with an electron distribution very close to 
the one derived from previous recurrence formulas. The 
Fermi-Dirac distribution is obtained in the limit of large 
S-values, with temperatures and Fermi levels that de- 
pend on the energy given initially to the system. Our 
computer program may handle single-electron level dis- 
tributions that could be difficult to analyze theoretically 
(for example, two bands of states). 

When the system is in thermal contact with a heat 
bath, electrons at level k are ascribed a probability per 
unit time, p, of being promoted to level fc-l- 1 and a proba- 
bility, qp, of being demoted to level k—1, provided these 
levels are empty. Strictly speaking, these prescriptions 
rest on an Einstein-type model of solids that supposes 
that the atoms are vibrating at frequency tOphonon = ^/fi, 
where e denotes as before the electronic level spacing. 
But the details of the thermalization model turn out to be 
rather unimportant. The computer program enabled us 
to reproduce the theoretical results with great accuracy. 
For example, when B = 100, e = ImeV, and T =100K, 
the numerical distribution fits the Fermi-Dirac distribu- 
tion with a discrepancy not exceeding 0.2%. 

When the electron gas is enclosed in a single-mode 
cavity, the probability that the cavity contains m light 
quanta is proportional to W{r — md), where W(r) de- 
notes the statistical weight of the electron gas for an 
added energy r, see Appendix ^ of If, initially, 

only the highest levels are occupied, we obtain ex- 
actly, from the recurrence relation satisfied by p{P', N, r), 
variance(TO)/average(m) — {B + l)/6. 
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APPENDIX B: OCCUPANCY AND 
FLUORESCENCE FOR TWO-SPIN-STATES 
ELECTRONS 

In the present Appendix we restrict ourselves to energy 
levels Sk = k, with k =^ 1,2..., the origin of the energy 
being set at A: = 0. Electrons are allowed to change spin 
in the course of time. We first consider an isolated system 
with constant numbers of spin-up and spin-down elec- 
trons (iV+ and N^, respectively), electrons of different 
spins being allowed to exchange energy. Next, spin- flip 
is allowed. Averaging, with the Boltzmann factor as a 
weight, provides occupancies for the case where the sys- 
tems are in contact with a heat bath. Occupancies in 
these various situations are illustrated in Fig. |l|, again by 
comparison with the FD distribution. 

Consider first an isolated system with iV+ spin-up elec- 
trons and N~ spin-down electrons, and suppose that the 
two sub-systems may exchange energy but that spin flip is 
not allowed. Setting N+ + N- = 2N and N+-N- = 2n, 
the system least energy is 



Un = 



N+{N+ + 1) N-{N- + 1) 

2 ^ 2 

N{N+l) + n'^ 



(Bl) 



If r > denotes the energy added to the system on 



top of N{N + 1), the remaining energy s 



splits 



into ri in sub-system 1 and r2 = s — ri in sub-system 2. 

We have shown in Appendix ^ that, for the case 
presently considered and in the limit — > oo, the num- 
ber of microstates for single-spin electrons is the num- 
ber of partitions p{r) of the excess energy r. The num- 
ber of microstates relating to one particular splitting 
of s is therefore p{ri)p{r2)- Accordingly, occupancies 
are obtained by averaging single-spin state occupancies 
shifted by ±n, with a probability law proportional to 
p{ri)p{s — ri) with ri running from to s. 

ne{r; k) = 

Z^r, ("i(r-i; fc - n)p{s - ri) + m{s - ri; fc + n)p{ri)) 



(B2) 



where m(r; k) is given by Eq. ( [aTiI ). Figure |i|c compares 
ne(r; fc) in Eq. ( p^ ) to npD in Eqs. (^), (||) for the case 
where n = {N^ — or r — s) and various values of 
the added energy. 

When spin flip is allowed, + N~ = 2N remains 
fixed, but n = {N^ — N~)/2 may take any value that 
does not exceed y/r, where r denotes as before the en- 
ergy added on top of N{N -1-1). If r = 6 for example, 
five values of n are permitted, namely n = 0, n = ±1 and 
n = ±2. It thus s uffi ces to sum the numerator and de- 
nominator of Eq. ( p3S| ) over permissible values of n. The 
occupancy reads 



ns{r; k) = 

J2ri k - n)p{s - ri) + m,{s - ri; fc + n)p{ri)) 



EnEri p{.s - ri)p{ri) 



(B3) 



where the sum over ri is from to s = r — n^, and m{r; fc) 
is given Eq. (All). The two terms in the numerator give 
equal contributions. Figure ^ compares ns[r\ fc) as given 
in Eq. (B3) with the FD distribution for various values 



of the added energy. 

When the system is in contact with a heat bath at 
temperature reciprocal /3, r fluctuates with a probability 
law where, as before, q = exp(— /?). Accordingly, occu- 
pancies are obtained by multiplying the numerator and 
denominator of the previous expression in Eq. ( |B3| ) by 
q^ and summing over r from to oo. The denominator 
gives the sum over states 



\ r / n 



(B4) 



which may be written as an infinite product |12|. The 
average added energy reads 



q dZ* 

E 

i=i,2.. 



2j 



1 



l3 - , 



(B5) 



The expression for the occupancy may be reduced to 
a double sum, which coincides with the one given in Eq. 
(^), shifted by ±n and weighted by . This final re- 
sult is intuitive since unbalancing between N'^ and 
increments the energy by n^. We have 

/ M En ("(g; k-n)+ n{q; fc + n)) g"' 
nc{q;k)^ — (B6) 

where the sums over n run from — oo to -l-oo, and n{q] fc) 
is given in Eq. (A_13). The two terms in the numerator 
give the same contributions. The occupancy nc{q; fc) may 
be expressed as a func tion of the average energy and k 
with the help of Eq. (p5|). Comparison with the FD 
distribution is in Fig. |l|e. 

Similarly the fluorescence is obtained by shifting fc by 
n with a weight factor g" . 
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FIG. 1. Ratio of exact and Fermi-Dirac (ufd) occupancies plotted as functions of ufd at equal values of the average added 
energy r. The isolated single-spin-state electron occupancy in a) was multiplied by 2 to account for the two spin states. In a) 
and b) spin flip is not allowed. In c) electrons of different spins may exchange energy but spin flip is not allowed. Cases d) and 
e) are the same as cases b) and c) respectively, except that spin flip is allowed, f) collects previous results for an average added 
energy r — 600. 
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FIG. 2. This figure is similar to Fig. |l|, but the comparison between canonical and grand-canonical ensembles is made at 
equal temperatures rather than at equal average energies. The energy level spacing is supposed to be 1 meV, typical of 1 /im 
long quantum wires. The parameter is the electron gas absolute temperature T. a) Spin flip is not allowed, b) Spin flip is 
allowed. 




10 20 30 40 10 20 30 40 

d d 

FIG. 3. Ratio of spontaneously-emitted light power (fluorescence) from 1 ^m-long quantum wires in contact with diamond 
(canonical ensemble) and copper (grand-canonical ensemble), respectively, as a function of d = huj/t. The parameter is the 
electron gas absolute temperature T. a) Spin flip is not allowed, b) Spin flip is allowed. 
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